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Abstract: The mainstay therapy against leishmaniasis is still pentavalent antimonial drugs; however, the rate of antimony 
resistance is increasing in endemic regions such as Iran. Understanding the molecular basis of resistance to antimonials 
could be helpful to improve treatment strategies. This study aimed to recognize genes involved in antimony resistance of 
Leishmania tropica field isolates. Sensitive and resistant L. tropica parasites were isolated from anthroponotic cutaneous 
leishmaniasis patients and drug susceptibility of parasites to meglumine antimoniate (Glucantime*) was confirmed using 
in vitro assay. Then, complementary DNA-amplified fragment length polymorphism (cDNA-AFLP) and real-time reverse 
transcriptase-PCR (RT-PCR) approaches were utilized on mRNAs from resistant and sensitive L. tropica isolates. We iden- 
tified 2 known genes, ubiquitin implicated in protein degradation and amino acid permease (AAP3) involved in arginine 
uptake. Also, we identified 1 gene encoding hypothetical protein. Real-time RT-PCR revealed a significant upregulation of 
ubiquitin (2.54-fold), and AAP3 (2.86-fold) (P<0.05) in a resistant isolate compared to a sensitive one. Our results suggest 
that overexpression of ubiquitin and AAP3 could potentially implicated in natural antimony resistance. 
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INTRODUCTION 

Leishmaniasis is one of the tropical diseases caused by infec- 
tion with the dimorphic protozoan parasites Leishmania. The 
clinical manifestations vary from self-healing cutaneous lesion 
to fatal visceral disease [1]. Cutaneous leishmaniasis (CL) is 
the most common form of the disease which has been report- 
ed from 77 countries with incidence of 1-1.5 million annual 
new cases [2], CL is an important public health problem in 
different parts of Iran and 13-25% of reported CL in some parts 
of Iran are anthroponotic cutaneous leishmaniasis (ACL) caus- 
ed by Leishmania tropica [3,4]- Chemotherapy still remains the 
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sole choice for treatment and antimonial compounds such as 
meglumine antimoniate (Glucantime®) have been in use for 
more than 70 years [5]. Unfortunately, in recent years antimo- 
nial therapy has been imperiled due to an increasing rate of 
treatment failure attributed to antimony resistant parasites in 
some endemic foci, especially in India where 37-70% of vis- 
ceral leishmaniasis patients failed to respond after antimonial 
therapy [6]. Failure of antimonial treatment was also reported 
in 10-12% of ACL patients in Iran which was shown to be due 
to drug resistance [7]. Most of the knowledge regarding anti- 
mony resistance process has derived from laboratory generated 
parasites suggesting presence of different mechanisms such as 
decreased capacity of converting SbV to Sblll, reduced drug 
uptake, and elevated drug efflux [8]. In laboratory generated 
parasites, elevated levels of trypanothione [9] and downregula- 
tion of aquaglyceroporin (AQP1) which uptake trivalent met- 
alloid [10] have been demonstrated. Additionally, overexpres- 
sion of ATP-binding cassette (ABC) protein MRPA which se- 
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questered Sb-thiol complexes has been related to the drug re- 
sistance [11]. 

Recent studies on resistant field isolates have suggested that 
in addition to the conventional mechanisms other factors such 
as apoptosis might be implicated in natural antimony resis- 
tance [12,13]. Considering that, more insight into the mecha- 
nisms responsible for natural antimony resistance could im- 
prove cuiTent treatment strategies. 

Several techniques such as RT-PCR microarrays, and pro- 
teomics are cunently used to explore molecular markers of an- 
timony resistance [14,15]. cDNA-amplified fragment length 
polymorphism (cDNA-AFLP) is one of the most sensitive tech- 
nologies for detection of differentially expressed gene which 
does not require any prior information of gene sequences [16, 
17]. In this study, cDNA-AFLP was used to identify the poten- 
tial biomarkers of antimonial drug resistance in L. tropica field 
isolates. 

MATERIALS AND METHODS 

Clinical isolates and cell culture 

L. tropica field strains used in the present study included 
meglumine antimoniate sensitive (MHOM/IR/10/Mash-175) 
and resistant (MHOM/IR/10/Mash-827) isolates which were 
obtained from patients with ACL in Mashhad, north-east of 
Iran [7]. Unresponsiveness was defined as patients who re- 
gardless of receiving at least 3 full courses of Glucantime® still 
showed active lesions [7]. The Institutional Ethical Committee 
of the School of Public Health, Tehran University of Medical 
Sciences reviewed and approved the study protocol. 

The parasites within 6 subcultures after the isolation from 
the patients were used. Cryopreserved parasites were recovered 
from liquid nitrogen and cultured in RPMI 1640 medium (Gib- 
co/BRL, Rockville, Maryland, USA) supplemented with 10% 
fetal bovine serum (FBS; Gibco/BRL), 100 U/ml penicillin, 100 
ug/ml streptomycin (Gibco/BRL), in the absence of antimoni- 
al drug pressure, and incubated at 26°C. 

In vitro drug sensitivity assay 

The drug susceptibility of the field isolates in monocyte-ma- 
crophage mouse cell-line (J774A.1) was analyzed in triplicate 
as describe previously with slight modifications [7]. Briefly, 
J774A.1 (5 x 10 4 cells/well) were grown in RPMI with 10% FBS 
in Labtek 8-well chamber slides (Nunc, New York, USA) and 
incubated at 37°C for 24 hr to allow cell attachment, then the 



cells were infected with stationary phase promastigotes at a 
parasite to macrophage ratio of 10:1. After 4 hr incubation at 
37°C, cells were washed for 3 times and incubated for 72 hr in 
the presence of serial dilution of Glucantime® (Sanofi-Aventis, 
Paris, France). Pentavalent antimony concentrations used for 
the sensitive isolate were 2, 4, 6, 8, 10, and 12 ug/ml, and the 
for resistant isolate were 35, 40, 45, 50, 55, 60, 65, and 70 ug/ 
ml (the doses were obtained based on previous screening test). 
The fresh drug was added, and slides were incubated for an 
additional 72 hr, and then stained using Giemsa. The number 
of amastigotes in 100 randomly chosen macrophages was 
counted. The inhibitory concentration 50% (IC 5 o) was calcu- 
lated using the values for the number of amastigotes/macro- 
phage. ICso was defined as the effective dose of Glucantime® 
that reduces the survival of L. tropica by 50%. 

RNA extraction 

Total RNA extraction was carried out using Tripure isolation 
reagent (Roch, Mannheim, Germany) according to the manu- 
facture's instructions with minor modifications. Briefly, 1 x 10 8 
log-phase promastigotes of sensitive and resistant isolates were 
harvested and lysed using Tripure reagent. The lysate was mixed 
with chloroform and vortexed vigorously, and then centrifug- 
ed at 12,000 g for 15 min at 4°C. Afterwards, the upper color- 
less phase was precipitated by adding isopropanol and washed 
with 75% ethanol. Finally, the precipitated RNA was dried at 
room temperature and dissolved in RNase free water. The quan- 
tity and quality of RNA were analyzed using nanodrop (ND- 
1000, Thermo Scientific Fisher, Waltham, Massachusetts, USA) 
and gel electrophoresis, respectively. To avoid any genomic 
contamination, RNA was treated with DNase (Qiagen, Hilden, 
Germany) according to the manufacture's protocols. 

cDNA-AFLP 

The cDNA-AFLP assay was conducted as describe previously 
with some modifications [18]. Single strand cDNA (sscDNA) 
was synthesized using 10 ug RNA, 20 pmol/pl oligo-dT (Fer- 
mentas, Burlington, Canada), 20 pmol/pl random hexamer 
(Fermentas), 10 mM of dNTP mix (Fermentas) incubated at 
65 °C for 10 min, followed by addition of 20 U RNase inhibitor 
(Fermentas), 4 pL of 5 x reverse transcriptase (RT) buffer con- 
taining Tris-HCl (pH 8.3) (Fermentas), and 200 U RevertAid 
premium RT (Fermentas). It was then incubated at 25°C for 10 
min followed by incubation at 50°C for 60 min. The integrity 
of cDNA was checked using glyceraldehyde 3-phosphate dehy- 
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drogenase (GAPDH) primers as the housekeeping gene (Table 
1). The PCR condition was an initial denaturing step of 94°C 
for 5 min and 30 repetitions of denaturation at 94°C for 30 sec, 
annealing at 60°C for 30 sec, and extension at 72°C for 45 sec 
with a final extension of 72°C for 7 min. 

Double strand cDNA (dscDNA) synthesis was performed 
using 40 U DNA Polymerase I (Fermentas), 1.5 U Ribonucle- 
ase H (Roche, Mannheim, Germany), and 10 mM of dNTP 
mix at 16°C for 2 hr, followed by purification with High Pure 
PCR Clean up Micro Kit (Roche). Double restriction digestion 
was carried out with 5 U Mbol and 5 U EcoRl (Fermentas) on 5 
ug purified dscDNA at 37 °C for 2 hr. 

Digested dscDNA fragments were ligated to AFLP adaptors 
(Table 1), 8 ug AD Mbol, AD EcoRl, 4 ug ad Mbol, and ad 
EcoRl in a final volume of 60 pi. Briefly, it was conducted in 
the following conditions: 1 min at 55°C, decreasing to 10°C 
over 45 min (i.e., 1°C per 1 min) and then 6 U T4 DNA Ligase 
(Roche, Mannheim, Germany) was added and incubated at 
4°C overnight. After purification with High Pure PCR Clean up 
Micro Kit (Roche), pre-amplification with pre-amp primers 
(Table 1) was performed according to the following conditions: 
30 cycles of 94°C, 30 sec, 64°C, 30 sec, and 72°C, 45 sec, and fi- 
nal extension at 72°C for 7 min. 

The sensitive amplification was done using various dilutions 
of pre-amplified products to get the optimal conditions. PCR 
amplification was carried out as described above using sensi- 
tive primers containing adaptor sequences plus 1 nucleotide at 



3' end (Table 1). 

The products resulting from sensitive amplifications were 
resolved on 8% non-denaturating PAGE and stained using sil- 
ver nitrate. Finally the gels were wrapped and scanned, and 
then checked for the presence of differentially expressed tran- 
scription-derived fragments (TDFs). 

Isolation, cloning, and sequencing of TDFs 

Selected bands were extracted from the polyacrylamide gels. 
The eluted DNA was re-amplified with selective primers (Table 
1) from resultant profile. After checking PCR products on 1.5% 
agarose gel, they were cloned into a pGEM-T Vector System I 
(Promega, Fitchburg Wisconsin, USA). To identify the isolated 
TDFs, recombinant plasmid containing an unknown DNA 
was subjected to sequencing using universal primers (T7 pro- 
motor and SP6) (Bioneer, Seoul, South Korea). Homology 
searches were performed in non-redundant nucleic and pro- 
tein databases BLAST (http://www.ncbi.nlm.nih.gov/BLAST). 

Real-time RT-PCR 

In order to confirm the differences in expression patterns of 
TDF-derived genes between resistant and sensitive L. tropica 
isolates, real-time RT-PCR was conducted. Target gene primers 
were designed by Primer 3 software version 0.4.0 (http://fro- 
do.wi.mit.edu) according to the identified genes (Table 1). The 
cDNA was synthesized from 1 ug of total RNA using Quanti- 
Tect® Reverse Transcription Kit (Qiagen, Hilden, Germany) 



Table 1 . Sequences of adaptors and primer used in cDNA-AFLP and real-time RT-PCR 



Name 


Sequences 


Adaptors and primers 


AD EcoRl 
ad EcoRl 
AD Mbol 
ad Mbol 


5'-ACCGACGTCGACTATCCATGAAG-3' 
5'-AATTCTTCATGG-3' 
5'-CACTATCCAGACTCTCACCGCA-3' 
5'-GATCTGCGGTGA-3' 


Adaptors 


Pre EcoRl 
Pre Mbol 


5-ACCGACGTCGACTATCCATGAAGAATTC-3' 
5-CACTATCCAGACTCTCACCGCAGATC-3' 


Pre-amplification primers 


51 EcoRl 

52 EcoRl 

53 EcoRl 

54 EcoRl 

51 Mbol 

52 Mbol 

53 Mbol 

54 Mbol 


5-ACCGACGTCGACTATCCATGAAGAATTCC-3' 

5'-ACCGACGTCGACTATCCATGAAGAATTCG-3' 

5'-ACCGACGTCGACTATCCATGAAGAATTCA-3' 

5'-ACCGACGTCGACTATCCATGAAGAATTCT-3' 

5'-CACTATCCAGACTCTCACCGCAGATCC-3' 

5'-CACTATCCAGACTCTCACCGCAGATCG-3' 

5'-CACTATCCAGACTCTCACCGCAGATCA-3' 

5-CACTATCCAGACTCTCACCGCAGATCT-3' 


Sensitive primers 


GAPDH F 
GAPDH R 


5-GCATGTGCTGACAAAGGAGA-3' 
5'-GGTCGTACTCGGGATGATGT-3' 


Housekeeping gene primers 


Ubiquitin F 
Ubiquitin R 
AAP3F 
AAP3R 


5-ATCACGCGCATCTCTCCTAC-3' 
5'-GTGGTGGTGTGGAGTCGAG-3' 
5'-GGTGTCATCTCGAACTGCTTC-3' 
5'-CACGAGGTAGATGACAGACATGA-3' 


Target gene primers 
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following the manufacturer's instructions. RT-PCR was per- 
formed in uiplicate with 20 pi volumes using SYBR® Premix 
ExTaq™ II (Takara, Tokyo, Japan) in a CFX96™ Real-Time Sys- 
tem (C1000™ Thermal Cycler) (BioRad, Hercules, California, 
USA). The PCR condition was as follows: activation at 95°C for 
30 sec, amplification at 95°C for 5 sec, 60°C for 30 sec for 40 
cycles. The relative value of the expression level of each gene 
was calculated by comparing the cycle threshold (CTs) of the 
target genes with that of the housekeeping gene (GAPDH) us- 
ing the 2 Mct method [19]. 

Statistical analysis 

In vitro drug susceptibility assay was presented as IC 5 o which 
was determined using the linear regression analysis. Each ex- 
periment was conducted at least thrice, and the results are ex- 
pressed as the mean ± standard enor of the mean (SEM). The 
significance of differences was determined using the Student's 
t-test. The level of statistical significance was P< 0.05. 

RESULTS 

Drug susceptibility of Leishmania tropica field isolates 

In vitro drug sensitivity of L. tropica to Glucantime® was de- 
termined in J774-1 cells. The results showed that IC 5 o values of 
Sb (V) for resistant isolate was 50.0 ±3.5 ug/ml, whereas for 
sensitive isolate was 6.0± 1.5 ug/ml. The ICso of resistant iso- 
late was about 8-fold higher than the sensitive strain. 

The cDNA-AFLP analysis of clinical isolates 

To identify differentially-expressed transcripts putatively as- 
sociated with natural antimony resistance, a cDNA-AFLP ap- 
proach was utilized on resistant and sensitive L. tropica field 
isolates. Total RNA was extracted and the sscDNA was synthe- 
sized from each sample. The integrity of cDNA was confirmed 
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Fig. 1 . Detection of GAPDH as the housekeeping gene (1 54 bp) 
by reverse transcriptase-PCR (RT-PCR) on agarose gel (1 .5%). M, 
50 bp (base pair) molecular weight marker; S, sensitive; R, resis- 
tant; NC, negative control. 



by specific primers of the housekeeping gene (GAPDH) (Fig. 1). 

The cDNA AFLP results revealed differential expressions of 
TDFs on 8% non-denaturating PAGE by silver staining using 
different primer combinations (Fig. 2). Differentially-expressed 
bands were extracted, re-amplified, and cloned into a TA-don- 
ing system. After then, the recombinant plasmids were sequen- 
ced. The sequencing results were searched in non-redundant 
nucleic and protein databases to find homologous sequences 
and gene identity. 

Blast analysis revealed that 1 clone highly matched with ubi- 
quitine implicated in protein degradation and another clone 
highly homologous to amino acid permease (AAP3) involved 
in arginine uptake. Moreover, a clone was found to be similar 
to a gene which coded a hypothetical protein. The details of 
sequencing results of identified genes are presented in Table 2. 

Validation of cDNA-AFLP analysis by real-time RT-PCR 

In order to validate the cDNA-AFLP results, the mRNA ex- 
pressions of ubiquitin and AAP3 in resistant and sensitive L. 
tropica field isolates were analyzed using the real-time RT-PCR. 
Fig. 3 shows a significant upregulation of ubiquitin (2.54-fold) 
and AAP3 (2.86-fold) in the resistant isolate compared to the 
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Fig. 2. Expression patterns of TDFs extracted from cDNA- AFLP 
PAGE. Sensitive amplification of cDNA-AFLP on a PAGE from 3 
different primer combinations of S3 Mbol/S2 EcoRI, S3 Mbol/Si 
EcoRI, and S4 Mbol/SA Mbol. Arrows point to differentially ex- 
pressed TDFs which were isolated with code numbers of 1 , 2, 
and 3 (associated genes which derived from these TDFs men- 
tioned in Table 2). M, 50 bp molecular weight marker; S, sensitve; 
R, resistant. 



Kazemi-Rad et al. : Ubiquitin and amino acid permease genes in association with antimony resistance in L. tropica 41 7 



Table 2. Differentially expressed transcription-derived fragments identified via cDNA-AFLP 



Code No. 


Accession No. 


Annotation 


Length (bp) 


E-value 


1 


XM-001 685021.1 


L. major amino acid permease (MP3) gene 


220 


1e-69 


2 


XM-003721 589.1 


L. major hypothetical protein gene 


215 


2e-34 


3 


Z1 4232.1 


L. major ubiquitin gene 


147 


5e-38 


3.5 r 


AAP3 


3 r 


Ubiquitin 
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Fig. 3. The relative expression levels of ubiquitin and amino acid permease (AAP3) genes in sensitive and resistant L. tropica isolates us- 
ing real-time RT-PCR. The expression of GAPDH was used to normalize the data. Values given are the mean±SD of 3 different experi- 
ments (P<0.05). S, sensitive; R, resistant. 



sensitive one (P< 0.05). 

DISCUSSION 

Emergence of antimony drug resistance has eroded treat- 
ment and control of leishmaniasis in endemic regions such as 
Iran [6,7]. The molecular mechanisms operating in natural 
antimony resistance are not well known. Different approaches 
are cunently used to elucidate potential factors implicated in 
this phenomena. In the current study, we utilized cDNA-AFLP 
to identify genes potentially involved in drug resistance of L. 
tropica field isolates. We found 1 unknown gene and also 2 
known genes including ubiquitin and AAP3 genes. Further- 
more, differential expression of known genes was confirmed 
by real time RT-PCR. Our findings showed a significant upreg- 
ulation of ubiquitine and AAP3 mRNA expression in the resis- 
tant isolate compared to the sensitive one. 

Ubiquitin belongs to the class of heat-shock proteins which 
play important roles in various cellular processes such as deg- 
radation of defective proteins, endocytosis, transcription, DNA 
repair, and apoptosis [20]. One of the well understood func- 
tion of ubiquitin is protein degradation via ubiquitin-protea- 
some pathway (UPP) by which it prevents cells from toxic ac- 
cumulation of abnormal proteins [21]. Through the UPP, ubi- 
quitin attached to lysine residues of the target proteins leading 



to degradation of the ubiquitin-tagged protein by the 26S pro- 
teasome [22]. Several studies have revealed that the wild-type 
ubiquitin is cytoprotective and vital for cell surveillance [21,22]. 
Tsirigotis et al. [23] studied the sensitivity of mammalian cells 
expressing mutant ubiqutin to protein-damaging agents such 
as cadmium suggested that genetic alteration within ubiqutin 
coding region inhibits protein degradation and sensitizes cells 
to toxic effects of abenant proteins. Another study showed that 
in mutant ubiquitin-expressing cells the viability was decreased 
when the cells were exposed to heat shock and hydrogen per- 
oxide [24]. 

It was demonstrated that metalloid compounds such as ar- 
senic and cadmium induce heat shock and/or oxidative stress 
responses leading to production of aberrant oxidized proteins 
[25,26]. Additionally, these compounds inhibit ubiquitination 
and degradation of aberrant proteins which might lead to oxi- 
dative stress-induced apoptosis [27]. It was documented that 
the elevated expression of ubiqutin in organisms such as Sac- 
charomyces cerevisiae contributes to resistance to oxidative stress 
through degradation of oxidatively damaged proteins [28]. 
Moreover, another study demonstrated that overexpression of 
ubiquitin in yeast increased ubiquitin-dependent degradation 
of damaged proteins and could prevent heath shock toxicity 
and cell death [29]. Therefore, our results suggest that the over- 
expression of ubiquitin in antimony resistant L. tropica isolate 
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could contribute to degradation of oxidized proteins and de- 
fend cells from oxidative stress related to antimonial com- 
pounds. 

Also, our results identified the upregulation of AAP3 in re- 
sistant isolate in comparison to the sensitive one. To date sev- 
eral AAPs have been identified in Leishmania [30]. AAP3 is a 
highly specific arginine transporter that belongs to a large ami- 
no acid/auxin permease family (AAAP) which localizes to the 
surface membrane of Leishmania [31]. Arginine is known as a 
precursor of polyamine biosynthesis [32] which transferred 
into Leishmania from external sources through AAP3 [31]. This 
amino acid is hydrolyzed by arginase producing omitine which 
participate in the synthesis of polyamine and trypanothione 
[32]. The trypanothione is a major reduced thiol of Leishmania 
parasites and plays a crucial role in detoxification of antimoni- 
al drugs [9,33]. Furthermore, the association of increased try- 
panothione levels with antimony resistance has been observed 
in Leishmania isolates [9,34]. Haimeur et al. [9] reported that 
the upregulation of trypanothione in resistant isolates is corre- 
lated with overexpression of enzymes involved in the synthesis 
of polyamines (ornithine decarboxylase; ODC). Based on the 
above data, it is concluded that an elevated expression of AAP3 
in antimony resistant L. tropica isolate might contribute to en- 
hance the level of trypanothione and consequently detoxifica- 
tion of antimonial drugs. 

In conclusion, for the first time, we identified ubiquitin and 
AAP3 genes as potential antimony resistance candidate target 
genes in Iranian stfains of L. tropica through cDNA-AFLP ap- 
proach. The overexpression of aforementioned genes might be 
implicated in natural antimony resistance and could be help- 
ful to get more insight into the molecular etiology of antimo- 
ny resistance in Leishmania spp. However, further studies are 
required to establish the exact mechanisms by which these 
genes act in the resistance process. 
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